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We study theoretically a dilute gas of identical fermions interacting via a p-wave resonance. We 
show that, depending on the microscopic physics, there are two distinct regimes of p-wave resonant 
superfluids, which we term "weak" and "strong". Although expected naively to form a BCS-BEC 
superfluid, a strongly-resonant p-wave superfluid is in fact unstable towards the formation of a gas of 
fermionic triplets. We examine this instability and estimate the lifetime of the p-wave molecules due 
to the coUisional relaxation into triplets. We discuss consequences for the experimental achievement 
of p-wave superfluids in both weakly- and strongly-resonant regimes. 
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Recently there has been considerable interest in trying 
to create a p-wave resonant superfluid experimentally. 
The BCS and BEG regimes for such superfluids are not 
just different aspects of the same phase, as they are for 
the s-wave resonant superfluids, but rather are different 
phases. Thus the tuning from the BGS to BEG regime 
involves a phase transition (or sometimes a sequence of 
phase transitions) [l|, Q ■ Such a transition can even be 
topological in some cases [H, S 0, [E @ • If the superfluid 
is confined to two dimensions, the BGS phase will be 
topological and will support vortices with non-Abelian 
excitations P, 

In this paper we show that resonant p-wave superflu- 
ids must be classified as two distinct types, with weak 
or strong Feshbach resonances (to be defined precisely 
later). The existing mean field theory of p-wave super- 
fluids, worked out in P, [3, applies only to the case of 
weak Feshbach resonances. However, as we shall establish 
below, the p-wave resonance used in ongoing experiments 
on 0] is a strong resonance. It is therefore im- 

portant to determine the properties of strongly-resonant 
p-wave superfluids. 

The full theory of strong p-wave resonances is yet to be 
constructed. Here we investi gate an effect first noticed by 
Y. Gastin and collaborators [lO|: in the regime of strong 
p-wave resonances the fermions form triplet states with 
angular momentum (spin) 1. Superficially similar to Efi- 



mov states these triplets are quite unusual. They are 
very strongly bound, with a binding energy largely inde- 
pendent of detuning from the resonance, as long as the 
detuning is not too large (but dependent on the strength 
of the resonance). Gorrespondingly, their size is of the 
order of the closed channel bound molecular state, far 
smaller than the average interparticle separation. We 
find the critical value of the resonance's strength at which 
the triplets first appear, and calculate their binding en- 
ergy as a function of the resonance strength. 

Thus if a BEG condensate of strongly-resonant p-wave 
molecules is created, one of its main channels of decay will 



be by molecular inelastic collisions, with two molecules 
turning into one atom and one triplet. We estimate the 
molecular lifetime due to this process and compare this 
with experimental observations We discuss limita- 
tions on the achievement of p-wave superfluids in both 
weak and strong resonances, arising from this and other 
inelastic decay processes. 

The theory developed here can be used to investigate 
the true ground state of a strongly-resonant p-wave con- 
densate. This is likely to be a gas of fermionic spin 1 
triplets (or possibly of larger composite particles). 

We consider a p-wave resonantly coupled superfluid, 
whose Hamiltonian is given by [l|, 0, 0] 
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Here , a are the creation and annihilation operators 
of a spinless fermion (atom) with mass m, and 6^ 
are the creation and annihilation operators of a bosonic 
diatomic molecule of spin 1 (the 3D vector index /i rep- 
resents the projection of spin). This superfluid is con- 
trolled by four parameters. The first two are the detun- 
ing eo and the overall particle number N , an expectation 
value of the operator TV = J^p^l^P + '^J^p^.qKJ'm- 
is more convenient to work with the energy equivalent of 

N, the Fermi energy ep = {Qtt'^K^N/V)^^'^ /{2m). The 
other two are contained in the coupling constant 5(|p|). 
The physical origin of the dependence oi g on |p| lies 
in the fact that the molecules have finite size. This can 
be captured by choosing g to remain constant as long 
as IpI A (which we denote simply by g) and quickly 
drop to zero if |p| ^ A. Here Re ~ h/K is the physical 
(closed-channel) size of the molecules. The knowledge of 
exactly how g drops to zero at large momenta may be 
important. In this paper we adopt the "hard momentum 
cutoff" approach 5(|p|) = .g6(A — p) (0 is equal to 1 or 



depending on whether its argument is positive or neg- 
ative). We have studied other types of cutoff, and find 
these do not change the main conclusions of this paper. 

Two dimensionless parameters can be constructed out 
of g, ep, and A, namely 



FIG. 1: The diagrams whose sum gives the scattering ampli- 
tude between an atom and a molecule 
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Notice that in order to observe universal (short distance 
physics independent) behavior, the interparticle separa- 
tion (~ h/ y/mep) must be kept much bigger than Re, 
thus 7 <C C2 . 

Both of these parameters control the perturbative ex- 
pansion of ^ in powers of the coupling g. It is cus- 
tomary, when analyzing Eq. ([T]), to apply a mean field 
approximation whose validity is based on the smallness 
of g. Strictly speaking, both 7 and C2 must be small 
in order for the mean field approximation ernployed in 
the original publications investigating Eq. ((T|) l^] to be 
valid. 7 depends on the interparticle separation and can 
be made small simply by reducing the particle density. 
C2 however depends solely on the physics of the Fesh- 
bach resonance which led to Eq. ([T]); its value, which can 
be small or large, is fixed by the atomic type and Fes- 
hbach resonance involved, so it cannot be continuously 
controlled. 

One terms the superfluids with 7 ^ 1 as those with 
narrow Feshbach resonances, while the ones with 7 3> 1 
are the broad Feshbach resonance superfluids [11, 0, . 
Likewise, we will term the C2 3> 1 resonances as the 
strong p-wave Feshbach resonances, while those with 
C2 <C 1 are weak resonances. The p-wave resonances are 
typically narrow because, even if they are not, they can 
be made narrow by reducing the particle density. 

The narrow and weak p-wave resonances have been 
thoroughly investigated in prior publications. It is there- 
fore imperative to consider the narrow and strong reso- 
nances. The main idea behind the analysis is based on 
the fact that fluctuational corrections to the mean field 
come from two distinct regions in momentum space, p 
of order h/l where I is interparticle spacing, and p of 
order A. The former capture the many-body physics of 
Eq. IT]) and are small as long as 7 is small. The latter 
come from high momenta and energies at which no real 
particles propagate. Thus this contribution, controlled 
by C2 , is essentially few-body, equivalent to solving some 
few-body Schrodinger equation, which although difficult, 
is not an impossible task. 

The analysis carried out in Ref. 0] showed that the p- 
wave Feshbach resonance in ""^K used in Refs. [1, [§| was 
strong. The derivation relied on the scattering amplitude 
of two atoms calculated in Ref. P, 0] 
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Here v is the so-called effective volume, controlled by 
the physical detuning [related to eo by loq — (eo — 
mg^ / (dn^h^)) / {1 + C2)] and ko is a parameter similar 
to the effective range of s-wave scattering (having, how- 
ever, the dimensions of inverse length). If fco and A are 
known (numerically or experimentally), C2 can be found 
from (g]). 

We also remark that had we considered a one-channel 
model of identical fermions interacting via a short range 
p-wave potential, such as V{r) = Xdfi6^^\r)d^, we would 
have obtained Eq. ^ with kg ~ — A Q. In other words, 
such a model automatically describes strong resonances. 

We now turn our attention to the physical conse- 
quences of strong resonances. Its main consequence is the 
existence of a bound state of three atoms when C2 exceeds 
a certain threshold. To show this, we calculate the scat- 
tering amplitude of one atom and one molecule. This 
is given by a sequence of diagrams depicted on Fig. [T] 
These diagrams are identical to the ones studied in the 
context of the s-wave BCS-BEC crossover 1^ 14 1. 

Here the atoms propagate with the free propagator 
G{p, uj) = 1/ (^Lo — ^ + iO^ , while to find the molecu- 
lar propagator one needs to calculate its self-energy Q 
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(In these, and subsequent, expressions we set h — 1 for 
clarity.) Each loop in the diagrams on Fig. [T]is linearly 
divergent. It is this divergence, occurring at momenta 
p ^ A and controlled by C2, which we would like to cap- 
ture. To do so, we study the atom-dimer scattering prob- 
lem with the following kinematics; a boson of spin /i and 
4-niomentum (0, k + E3) scatters off a fermion with 4- 
momentum (0,0). The outgoing particles are a boson 
with spin 1/ and 4-momentum (q, go + k + E3) and a 
fermion with (— q, —qo)- Here k{ujo) is an implicit func- 
tion of the detuning such that the bosonic propagator 
D(q, go + k) has a pole as q, go ^0. i?3 < is the 
energy at which we are looking for a bound state. The 
scattering T-matrix has the following general form 

T^i^{p,Po) = Ti{p,pq)5^u + T2{p,PQ)py,Pulp^ (6) 
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and the scattering length Obf is related to Ti(0,0) (eval- 
uated at E3 = 0) as abf = ^Ti{0,0). 



J 



The integral equation for the T-matrix is derived anal- 
ogously to the s-wave problem 13|, , and is 



2 f d^q 

T^Ap,Po) = G{p,po + K + E3)p^p,g{\p\)g{\p\/2)-Ai —i-T^^{q,qo)D{(i,qo + K + E3)G{-(i,-qo) 

1 + C2 J [^TTr 



1 + C2 

xG{p + q, po 



go + K + E-i)ip + q/2)^{q + p/2),g{\p + q/2|).9(|q + p/2|) 
I 



(7) 



The factor I-I-C2 is the inverse residue of the bosonic prop- 
agator. Tpi/((f, qo) is analytic in the upper halfplane of go 
and thus we may integrate out go, setting go — g^/2m. 
To solve the integral equation we then let po —>■ —p^ /2m. 
For simphcity define Ti{p,—p'^/2) = Ti{p). Measuring 
momenta in units of the cutoff, energies in units of IS?' jm^ 
and the T-matrix itself in units of 1/(toA) we find the in- 
tegral equation 



67r 



2 C2 



1 + C2 p2 - K - 



-e(i-p)<52, 



-3C2 / q^dqD{q,-q^l2 + n^E^)a,,[j>,q)T,[(i) (8) 

The coefficients aji (p, g) are given by an integration over 
directions of q 
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FIG. 2: (a) Scattering length abf in units of fij tv and (b) 
binding energy Ez of the triplet in units of 1^ /m, both as 
functions of C2. Here, detuning has been set to zero. The 
large C2 limit is indicated. 
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The scattering length abf is found by solving Eq. ^ 
at i?3 = 0. The binding energy of the triplet corre- 
sponds to a pole in the T-matrix and thus to a solution 
of the homogeneous integral equation at a specific value 
of i?3. Fig. [2lja) shows how the scattering length is nega- 
tive for a weak Feshbach resonance, becoming more neg- 
ative and diverging at C2 w 3.3. This is the strength 
of the resonance at which the bound triplet appears, 
as illustrated in Fig. [D^b). As C2 ^ cxd the scattering 
length saturates at Obf ~ I.9/A and the binding energy 
at E^ « — O.llA^/m. The existence of a finite C2 ^ 00 
limit can also be seen by observing that each of the di- 
agrams depicted on Fig. [T] has a finite C2 ^ 00 limit, ft 
should be noted that whereas the presence of the bound 
state and the general features of the scattering length and 
binding energy do not depend on the method of cutoff, 
the exact numerical values will in general depend on the 
method chosen. 

The existence of the bound trimcr state for large 
C2 raises the possibility of an inelastic decay channel 
in which two dimers collide to leave a trimer and an 



unbound atom (with large relative velocity). (Hence- 
forth we use the term "dimer" to refer to a molecule 
of two atoms, to distinguish this clearly from a triplet, 
or trimer.) In a non-degenerate gas of dimers, these in- 
elastic losses will cause the density of dimers rid to decay 
as 
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with add — 2|j(fcitTin(fci)) where the average is over 
the relative momenta of the incident dimers, fc^, and 
o'in(fci) = / \fin\^dQ.^ with /in thc iuclastic scattering 
amplitude into a final momentum kf. Arguments similar 
to the ones presented above for dimer-atom scattering 
show that l/inP ^ Rl- Thus, for large C2, such that the 
trimer binding energy is —E3 ~ / (mR^) and is large 
compared to the incident kinetic energy, one finds 
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It is instructive to compare this result with thc inelas- 
tic decay constants into deep bound states for s-wave 
dimers, formed from (two-component) fermions or from 
bosons with s-wave scattering length a. Close to thc s- 
wave resonance, a ^ Re, the decay constant (|lip is much 
smaller than that expected for bosons, add*^"*^"" fia/m, 
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but is larger than that for s-wave dimers of fermions, 

„ s— fermion 
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15| . The suppressed decay of 
s-wave dimers of fermions is explained in Ref . [l5| as an 
effect of the Pauli principle, reducing the probability to 
find three atoms within a lengthscale i?e. In a p-wave 
dimer the two atoms have a probability of order unity 
to be inside the centrifugal barrier, at a separation of 
order Re. Taking this feature of the p-wave dimers into 
account, simple estimates lead to add ~ TiRe/m for de- 
cay into trimers, consistent with the result from the 
T-matrix calculation. In addition to this channel, there 
are inelastic channels ~ active for both weak and strong 
resonances - involving decay into deep dimer states. Ap- 
plying the same simple estimates, one finds that the in- 
elastic decay constants for dimer-dimer and dimer-atom 
scattering are also add ~ Q^da ^ TiRe/m. 

In recent experimental work a gas of p-wave Fesh- 
bach dimers was created in ^°K. Unfortunately the life- 
time of the dimers was observed to be quite short, about 
2 ms. While ""^K can suffer losses through dipolar relax- 
ation (an effect expected to be absent for p-wave res- 
onances in other fermionic systems, for example ^Li), 
Ref. ^] found that the lifetime was shorter than that 
predicted for dipolar relaxation alone. Additional losses 
could arise from inelastic collisions of the dimers. This 
mechanism would imply a density dependence of the de- 
cay rate; this dependence has not, as yet, been estab- 
lished experimentally. 

Within the above considerations, we expect the decay 
rate of dimers via relaxation into deep trimers or dimers 
under inelastic collisions (with other dimers or with un- 
bound atoms) to be of order ri,, ~ ^^n, where n is the 

y 111 1 

density of atoms or dimers with which a given dimer can 
collide. Taking n ~ 7 x lO^^cm"'^ (the atomic density 
in the experiments of Ref.[^) we find Fin ~ lOHz. This 
estimate is more than one order of magnitude smaller 
than the additional decay rate required to account for 
the observations of Ref. [9|. However, we note that the 
prefactor to the estimate is uncertain. In view of this 
uncertainty, and in view of the lack of clear evidence of a 
density dependence in the experiment, it remains an open 
issue whether the dimer lifetime in Ref. 9] is limited by 
inelastic collisions. 

Our analysis has important consequences for possibil- 
ities to achieve superfluid phases close to a p-wave reso- 
nance. On the BEG side of the resonance, our calcula- 
tions show that the elastic dimer-dimer scattering ampli- 
tude is /oi Re- Consequently, the elastic scattering rate 
is of order Fei '-^ ^^^ndihRe) , which is typically much 
smaller than the inelastic decay rate, Fjn ~ ^^rid. (For a 
BEG of dimers, ki is small compared to the inverse parti- 
cle spacing, l/l, so kiRe ^ Re/l <^ 1-) It is therefore un- 
likely that a BEG of dimers can undergo sufficient elastic 
collisions to thermalize before inelastic losses deplete the 
gas. On the other hand, on the BCS side of the resonance, 
thermalization can proceed at a much faster rate, and 



will be limited by the rate of hybridisation of the dimers 
with the unbound atoms (this is the rate at which pairs of 
atoms can exchange their relative momentum). Assum- 
ing the densities of dimers and atoms to be comparable, 
n-d ^ n-ii = n, one finds that the hybridisation rate, as 
set by the width of the resonance ([3]), is Fhyb ^ j^^- 
Using ([2]), we find that Fhyb/rin ~ Thus, provided 

the resonance is not very weak (c2 very small), the rate 
of hybridisation is parametrically the same as Fjn, and 
the system may thermalize before inelastic losses deplete 
the gas. Thus our results show that it is on the BCS side 
of a strong resonance that one has the best opportunity 
to attain a thermalized p-wave superfluid phase. Finally, 
we note that the limitations we have described in this 
paragraph, arising from decay into deep bound states, 
could be eliminated in an "optical Feshbach" scheme in 
which the particles are coupled to a deep closed-channel 
molecule. In this case, it is important that the resonance 
should be sufficiently weak in order also to eliminate in- 
elastic decay processes into the triplet states that always 
exist for strong resonances. 
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